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SUMMARY 


The  diffraction  of  a  plane  electromagnetic  wave  by  a  thin  dielectric  wedge  is  studied 
using  a  set  of  approximate  boundary  conditions  together  with  the  Wiener-Hopf  technique. 
An  integral  expression  for  the  scattered  Field  is  derived  and  the  diffracted  far  field 
obtained  by  application  of  the  method  of  steepest  descents.  The  associated  GTD 
diffraction  coefficients  are  exact  to  order  O(A^),  where  2 A  is  the  (small)  wedge  angle. 
The  results  are  useful  for  wedges  for  which  i  A  «  1,  where  n  is  the  refractive 

index  of  the  dielectric  relative  to  that  of  the  surrounding  medium.  v 
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1  INTRODUCTION 

At  present  there  is  great  interest  in  the  electromagnetic  scattering  properties  of 
various  dielectric  and  radar  absorbent  materials  (RAM).  This  information  is  required,  for 
example,  for  the  calculation  of  the  radar  return  from  dielectric  and  RAM  coated  bodies. 
An  important  aspect  here  is  the  study  of  diffraction  by  the  edges  of  such  bodies.  In  the 

geometrical  theory  of  diffraction  (GTD)^  a  wave  diffracted  by  an  edge  is  characterised  by 

a  diffraction  coefficient.  For  high  frequencies  the  most  significant  contribution  to  the  edge 
diffraction  coefficient  may  be  determined  from  the  solution  to  the  problem  of  the 
scattering  of  a  plane  w'ave  by  an  infinite  wedged  In  the  absence  of  a  general  solution  to 
the  problem  of  plane  wave  diffraction  by  an  infinite  pentrable  wedge  there  have  been 

several  approximate  formulations.  One  such  formulation  is  to  consider  the  problem  of 
diffraction  by  a  thin  dielectric  w edge.  Kaminetzky  and  Keller^  have  derived  a  set  of 
approximate  boundary  conditions  under  the  assumption  of  a  small  wedge  angle.  The 

leading  contributions  to  the  diffraction  coefficients  were  then  deduced.  Analogous 
boundary  conditions  have  been  used  by  Leppingtom\  in  a  study  of  travelling  waves  in  thin 

dielectric  slabs,  and  by  Anderson4  and  Chakrabarti^,  in  their  w'ork  on  diffraction  b\  a 

dielectric  half-plane 

In  the  present  work,  approximate  boundary  conditions  are  used  in  conjunction  with 
the  Wiener-Hopf  technique  to  derive  an  approximate  solution  for  the  scattered  field 

external  to  a  thin  dielectric  wedge.  An  expression  for  the  diffracted  far  field  is  then 
obtained  by  the  steepest  descents  technique.  The  arsociated  diffraction  coefficients  are 

exact  to  order  O(A^),  where  2A  is  the  (small)  wedge  angle.  The  derivation  of  these 

coefficients  extends  the  results  contained  in  reference  2. 

In  the  following  section  we  formulate  the  problem  of  the  scattering  of  a  plane  wave 

incident  normally  upon  an  infinite  dielectric  wedge.  An  outline  of  the  derivation  of 

approximate  boundary  conditions  applicable  to  a  wedge  of  small  half-angle  is  then  given 
The  next  tw'o  sections  deal  with  the  details  of  the  Wiener-Hopf  method  of  solution  and 
with  the  extraction  of  expressions  for  the  diffraction  coefficients.  Finally  the  limitations  of 
the  analysis  are  discussed  and  the  results  are  compared  with  those  of  other  authors  and 

illustrated  by  numerical  examples. 
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FORMULATION  OF  THE  PROBLEM 


Consider  the  electromagnetic  scattering  of  a  plane  wave  incident  normally  upon  an 
infinite  homogeneous  dielectric  wedge.  The  wedge,  of  half-angle  A,  has  permittivity  (\ 
and  permeability  fj.0%  and  is  surrounded  by  another  dielectric  medium,  the  latter  having 
permittivity  (0  and  the  same  permeability  fiQ.  At  different  times  rectangular  cartesian 
coordinates  (x,y,z)  and  circular  cylindrical  coordinates  (r,0,z)  are  used,  with  x  =  rcosfl, 
y  =  rsin0  and  -x  ^  6  £  %.  The  z-axis  is  chosen  to  lie  along  the  edge  of  the  wedge. 
The  interfaces  between  the  two  media  are  at  0  =  ±(x-A),  as  shown  in  Figure  1.  The 
incident  plane  wave  approaches  the  wedge  from  the  0  =  0O  direction  (we  assume  0  ^  0O 
<x-A)  and  is  assumed  time-harmonic.  All  fields  will  then  exhibit  the  same  time 

dependence,  chosen  to  be  e-10^,  which  will  be  suppressed  throughout. 

Without  loss  of  generality  we  confine  the  analysis  to  the  cases  where  the  plane  wave 
is  either  H-polarised  or  E-polarised.  For  an  H-polarised  incident  wave  the  total  H-field 
everywhere  has  a  non-zero  z-component  only.  For  example,  the  external  solution  is  of 
the  form 


H  «  (0,  0,  Hzo(r,0)) 

-iu,'fQE  -  Z  x  H 


with 


where* 


H 

20 


(r ,  fl) 


0 


and 


o 


2 

(J)  u  ( 

o  o 


(la) 

(lb) 


(2) 


Analogous  expressions  describe  the  electromagnetic  field  within  the  wedge,  where  we 
denote  the  z-component  of  the  H-field  by  Hzj(r,0)  and  the  wavenumber  k0  is  replaced 
by  k| ,  defined  by  k  J  =  cxj^^0f].  Continuity  of  the  tangential  components  of  the  electric 
and  magnetic  fields  at  the  wedge  surfaces  gives: 


Hzo(=  (.-A)) 

-  Hz](s  (.-A)) 

(3a) 

<'-A» 

-  '<AH21<!  <"*>>  ■ 

(3b) 

Here  the  r-dependence  of  the  Hzo  and  Hzj  functions  has  been  suppressed. 

Similarly,  for  an  E-polarised  incident  plane  wave  the  total  field  external  to  the 
w'edge  is  given  by 


E  -  (0,  0,  Ezo ( r , 0 ) )  (4a) 

iaj^0  H  -  Z  x  E  .  (4b) 

where  Ez0  satisfies  the  wave  equation 

_  K +  -il  yr'«  -  0  <5> 

‘Throughout  the  paper,  3X  denotes  a  partial  derivative  with  respect  to  the  variable  x. 
The  extension  to  higher  order  derivatives  is  obvious. 


2 


In  this  case  the  boundary  conditions  take  the  form 


Ezo(i  (*-A) )  -  Ezl ( t  (»-A)) 

3f)Ezo(i  (*-A))  -  3eEzl(i  (tr-A)) 


(6a) 

(6b) 


where  Ezj  is  the  z -component  of  the  electric  field  inside  the  wedge. 

To  avoid  duplication  of  the  subsequent  analysis  it  is  convenient  to  consider  scalar 
fields  anc*  defined  outside  and  within  the  wedge  respectively,  satisfying  the  wave 
equations 


R  +  >3]  -  0  (7a) 

and 

R +  kil  ■  0  (7b) 

in  their  respective  domains  and  the  boundary  conditions 

O-A))  “  a  v^l  ( ~  (»-A))  (8a) 

dflSot*  (r~A))  "  b  (»-A))  .  (8b) 


With  an  H-poIari$ed  incident  wave  s  Hz  and  a  =  1,  b  =  <  ^  e  j.  For  the 
E-polarisation  case  v~  =  Ez  and  a  =  b  =  1.  The  incident  wave  is  chosen  to  have  unit 
amplitude.  Therefore  we  write 

i  nc 

v'  -  exp  [  -  i  k  (xcosP  +  vsin^  )]  (9) 

o  o  o  o 

We  are  primarily  concerned  with  obtaining  a  solution  external  to  the  wedge.  In  this 
region  we  write  =  v onc  +  v,  so  that  v  represents  the  scattered  field  and  satisfies 

R  +  o)  v(r,p)  *  0  (i0) 

At  present  a  general  solution  to  this  problem  is  not  available.  However,  if  the  wedge 

half-angle  is  assumed  small  (A  <<  1),  we  can  turn  to  a  perturbation  analysis 

3.  DERIVATION  OF  APPROXIMATE  BOUNDARY  CONDITIONS 

In  this  section  we  derive  a  set  of  approximate  boundary  conditions  (abc’s)  applicable 
to  a  study  of  scattering  by  a  thin  dielectric  wedge.  These  conditions  take  the  form  of 
expressions  for  the  discontinuities  in  the  external  field  and  its  normal  derivative  across  the 
half-plane  formed  by  the  wedge  faces  in  the  limit  A  ->  0. 

With  A  «  1,  Taylor's  theorem  may  be  used  inside  the  wedge  to  give 
v^(ir-A)  -  ^(-tt+A)  -  -A(8^j(*-A)  +  (~*+A) )  +  O(A^)  (11a) 

(*‘ A)  (-*+A)  -  -A  (3*^  (  -tt  +  A)  )+0(  A3  )  (11b) 
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Next,  the  (exact)  boundary  conditions  (S)  are  used  to  eliminate  all  -dependence  from 
these  relations.  This  gives 


b[v"  (*-A)*v  )]  -  -Aa[d  (  t-A)+8  f  ( -*  +  A)  ]  +0(  ) 

o  o  r  o  #7  o 

2  2  2  2 

arc*  v?  (ir-A)-8_*?  ( -  if  +  A )  ]  -  Ab  •  r  8  f  (*-A)+r  8  o  (-*+A) 

‘  o  fi  o  rr  o  rr  o 

+  r  8  ^  ( ?r  -  A)  +  r  8  vr>  (-ir+A) 


(12a) 


(12b) 


+  k^  r^vrQ(T-A)+  kj  (-x+A)  ]  +  O(A^) 

where  we  have  used  the  wave  equation  (7b)  to  rewrite  (lib)  in  a  form  suitable  for  the 
application  of  exact  boundary  conditions. 

It  is  unlikely  that  a  solution  can  be  found  for  v-0  subject  to  the  abc's  (12)  on  the 
wedge  faces.  However,  has  been  completely  eliminated  and  from  now  on  we  need 

only  consider  the  region  external  to  the  wedge.  Since  the  wedge  angle  is  small  this 
region  is  w'ell  approximated  by  all  of  three-dimensional  space  except  for  the  half-plane 
y=0,  x<0,  across  which  discontinuities  in  an^  e*ist-  Mathematically,  this  is 

justified  by  further  application  of  Taylor’s  theorem  provided  that  the  domain  of  is 
extended.  Using  also  (7a)  the  abc's  become 


b:.-0: 

-  -A(a-b)  (3^*0(*)  +  ) 

(13a ) 

and 

a;Vo] 

-  Ablkrto)  r2<‘-o(I>  *  >■-„<-'» 

♦  »<a-b)(ejf,.0<.)  . 

(13b) 

where 

(v*o!  E  v-0(T>  - 

(14a'1 

(3p*-o;  *  dfis0(z)  -  dfiv-p(-z) 

(14bN. 

and  the  omitted  terms  are  at  least  of  order  O(A^). 

In  terms  of  v  the  abc's  take  the  following  form  (in  cartesian  coordinates) 
b(v(x,0+)  -  v'(x.O-))  +  A(a-b)  x  (8  v  ( x .  0+)  +  d^v  (x,0-)) 

-ik  x  cosP 

2A(a-b)  ik  x  sin/>  e  °  °  ;  x  <  0  ( 1 5  a  > 

o  o 
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a(d  v(x,0+)  -  3  v'(x.O-))  -  Ab(k?-kh  x(v^(x ,  0+)  +  tf<x,0-)) 

y  y 


-  A(a-b)  [x3^  v'(x,0+)  -  3x^(x,0+)  +  x3^  tf(x,0-)  -  3^(x,0-)j  - 


_  *  -ik  xcosfl 

2Ab[krkJ  x  e 


.  -ik  xcosd 

+  2A(a-b)  ik  ( cosd  +  ik  xs i n  6  )  e  °  °;x<0.  (15b) 

o  o  o  o 


4  METHOD  OF  SOLUTION 

The  problem  has  been  reduced  to  one  of  solving  the  wave  equation  (10)  for  the 
scattered  Field  ^  subject  to  the  abc*s  (15)  and  the  radiation  condition,  together  with 
appropriate  edge  conditions.  In  addition  the  following  conditions  hold  on  the  y=0  plane 
for  x>0: 

^(x,0+)  =  v'(x.0“)  I  x  >  0  (16a) 

3yv(x,0+)  =  3yv(x,0-)  ;  x  >  0  .  (16b) 

The  form  of  the  relations  (15)  and  (16)  suggests  a  Wiener-Hopf  type  approach  for 

solution. 

Define  the  Fourier  transform  of  v(x.y)  (with  respect  to  the  variable  x)  by 

CO 

♦(s.y)  -  |  v'(x.y)  e*SX  dx  -  ^(s.y)  +  ’t_(s,y)  ,  (17a) 


where 


♦+(s,y)  -  [  v’(x,y)  e,sx  dx  ,  (17b) 

0 

0 

♦_(s,y)  -  |  ^(x.y)  e,sx  dx  .  (17c) 


Let  kQ  have  a  small  positive  imaginary  component  (kQ  =  k£  +  ik^,  0  <  k^  <<  k£)  so 
that  the  transform  ♦  exists  and  is  analytic  in  the  strip  -klQ  <  Im  s  <  k|j  cos^Q. ^  (The 
assumption  of  a  slightly  lossy  medium  means  that  Sommerfeld’s  radiation  conditions  are 
equivalent  to  the  requirement  that  the  scattered  field  vanishes  at  infinity*^.)  and 

are  analytic  in  the  overlapping  half-planes  Im  s  >  -kj,  and  Im  s  <  k^  cosfl0  respectively. 

Taking  the  Fourier  transform  of  (10)  and  solving  the  resulting  differential  equation 

gives 


s 


<Ks.y) 


B,(s)  e">(S):' 

;  y  >  0 

Bj(.)  »-'<s))' 

;  y  <  0 

(18) 


Vk-ith 


^(s)  "  I*2  -  ko) 


The  s-plane  is  cut  by  straight  lines  from  s  =  kQ  to  infinity  in  the  upper  half-plane  and 
f~om  s  =  -kQ  to  infinity  in  the  lower  half-plane.  The  branch  of  the  square  root  for 
which  7(0)  =  -ik0  is  chosen.  Upon  transforming  equations  (16)  for  x  >  0  w^e  find 


*+(s,0+)  -  *+(s,0-)  e  *+(s,0) 

dy*+(s’0+)  -  3yi+(s,0-)  e  3y*+(S,0) 


(19a) 

(19b; 


Similarly,  transforming  equations  (15)  for  x  <  0,  and  using  (19),  gives 

bii<s,0+)  -  ^(s.O-))  -  iA(a-b)f3^  ^(s,0+)  4  3^  ^(s,0-)]  4  bP  (s) 

l  s  v  s  v  J  4 


2A(a-b)  ik  sinfl 
o  o 

~  2 
( s - k  cos P  ) 
o  o 


(20a> 


(  2 

a  1 3  i(s,0+)  -  3  i(s,0-))  4  iAb  k  -k  |(3  i(s,0+)  +  0  ^(s,0-))  + 
y  y  l  1  o  J  s  s 

iA('a-b)f3^  i(s,0+)  -  s'k(s,0+1)  4  3^  ^(s,0-)  -  si(s,0-)l  4  aQ  (s) 

L  syy  syy  J  + 

2  A  T  b  f  k ^ 1  +  ( a  -  b )  k  fscos^  - k  )] 

'  L  I  o  J  o  00 


(2010 


where 


(s-k  cos^  )* 
o  o 


bP  (s)  -  2A(a-b)  i  a2  i  (s.O) 

+  SV  4 


( 2 1  a  l 


aQ+(s) 


2Ab  i f kf -k2 1  6  +  (s.O)  - 

11  0;  S  4 


2A(a-b)  i  (s  i+(s,0)  4  7  (s)  3s'k+(s,0)) 


(21b) 
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Applying  (18),  these  equations  may  be  rewritten  as 


f (s)  F_(s)  +  P+(s)  +  C(s)  -  0 

(22a) 

gCs)  C.(s)  +  Q+(s)  +  DCs)  -  0 

(22b) 

where 

F.(s) 

-  Bj (s)  -  B2(s)  -  ^ ( s , 0+ )  -  ¥_(s,0-) 

(23a) 

C_(s)  - 

-7  (s)  (Bj  (s)  +  B2(s))  -  3y'J'_(s,0+)  -  3y*_(s,0-) 

(23b) 

r<s)  -  1  ♦  1A  <abb)  -  >(.)  ££>] 

(24a) 

g(s> 


1  + 


i  Ab 


a  ■)  (s) 


C_(s) 


■)2(s) 


C_(s) 


-  iA^7(s)  §=$•  (24b) 


C(s)  -  -2  i  A 


(a-b) 


k  s  i  n  0 
o  o 


k  (s-k  cosfl  )  ^ 
o  o 


(25a) 


D(s)  -  - 


~A  [ bk^  -  ak^+  (a-b)  sk  cos£  ] 
2A  t  o  o  oJ 


(25b) 


(s-k  cos^  ) 
o  o 


and  the  prime  denotes  the  derivative  of  a  function  with  respect  to  its  argument.  (22)  are 
independent  Wiener-Hopf  equations.  Standard  techniques  may  be  used  to  solve  equations 
of  this  form^  for  the  unknown  functions  F_,  G_,  P+  and  0+  (regular  in  the  appropriate 
half-pianes)  provided  that  f,  g ,  C  and  D  are  known  functions  of  s,  regular  in  the  strip 
<  1m  s  <  cos^0.  In  the  current  application  C  and  D  have  been  determined,  but 
f  and  g  depend  on  the  unknown  functions  F_  and  G_. 

Suppose  for  the  moment  that  the  Wiener-Hopf  equations  may  be  solved  for  F_  and 
G_.  Then  the  scattered  field  C(x,y)  may  be  determined  from  the  Fourier  inversion 
formula 


v (x,y) 


oo+  i  c 


--XH-  i  c 


-yiyi-isx  fC,(  s) 
t-)(s) 


ds 


(26) 


where 

-k'<c<k1cos^ 
o  o  o 

Let  us  now  proceed  formally  with  the  method  of  solution  of  equations  (22)  proposed 
by  Jones^.  We  assume  factorisations  of  the  functions  f  and  g  of  the  form 

f( s)  -  f+(s)  F_(s)  (27a) 

g(s)  -  g+(s)  g_(s)  .  (27b) 
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Q  (s) 

K(s)  =  TTs)  +  T+(s)  *  ‘B-(s)  C-(s)  ’  T_(s)  •  (30b' 

These  equations  serve  to  define  J(s)  and  K(s)  in  the  strip  -k'0  <  Im  s  <  cos';0.  We 

note  that  the  second  part  of  each  equation  is  defined  and  regular  in  the  region  lm  s  > 
while  the  third  part  is  defined  and  regular  in  the  region  Im  s  <  k^  cos/0.  Since 
these  two  half-planes  overlap  analytic  continuation  may  be  used  to  deduce  that  the 
functions  J  and  K  are  defined  and  regular  in  the  whole  of  the  complex  s-planc.  If.  in 
addition,  J(s)  and  K(s)  tend  to  zero  (and  have  algebraic  behaviour)  as  s  tends  to  infinity 
in  all  directions  it  follow's  from  an  extended  form  of  Liouville‘s  theorem^  that  these 
functions  are  identically  zero.  Such  conditions  on  the  asymptotic  behaviour  of  J  and  K 
are  necessary  to  ensure  a  unique  solution  to  the  problem  A  unique  solution  will  be 
obtained  once  suitable  conditions  on  the  behaviour  of  v  at  the  edge  x  =  y  =  0  have  been 
imposed. 

To  specify  a  sufficient  set  of  edge  conditions  we  must  determine  the  behavioui  of 
the  quantities  in  the  second  and  third  parts  of  (30)  as  s  tends  to  infmitv  in  the 
appropriate  half-planes.  Let  us  suppose  that 
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V  ( x ,  y ) 

O(x^)  as  x  -*  0+  on  y  -  0 

(31a) 

v (*.y) 

0(x1')  as  x  ^  0-  on  y  -  Ot 

(31b) 

vv(x,y) 

^  0(\P)  as  x  -4  0+  on  y  -  0 

(31c) 

Vy(x.y) 

^  0(xa)  as  x  ->  0-  on  y  -  Oi 

( 3  1  d) 

It  follows  from  Abelian  theorems  on  Fourier  transforms  that 


F  (s^  ^  0(s  1  as  s  of  in  Im  s  <  k*cos^ 

o  o 


C  ^s'j  ^  0(s  °  *  )  as  s  -*  oo  i  n  I  m  s  <  k'cos/? 

o  o 


P  (s)  ^  0(s  ^  as  s  -*  w  in  Im  s  >  -k ^ 
+  o 


Q  (s)  ^  0(s  ^ )  as  s  -»  <*  in  Im  s  >  - k 
+  o 


From  (25),  (29a)  and  (29c)  is  is  clear  that 


S  (s)  -*  0  as  s  ■*  *  in  Im  s  <  k  cos 0 

o  o 


T  (s)  ->  0  as  s->tf  in  Im  s  <  k  cos*5 

oo 

It  remains  to  determine  the  asymptotic  behaviour  of  f*(s)  and  g+(s). 
(32a)  and  (32b)  we  find  that 


f(s)  l  +  i  Ac  j  sgn^s^  as  isi 
g(s)  ->  1  +  iAc2  s grids')  as  isi 


(32a) 

(32b) 

(32c) 

(32d) 

(33a) 

(33b) 

From  (24a),  (24b), 

(34a) 

(34b) 


where  the  constants  cj  (i  =  I  .2)  are  given  by 

r ( a-b) 


(35a) 


(rr  +  1>  ( a-b) 

c  -  - 

2  a 

and 

♦  1  for  Re  s  \  0 


(35b) 


sgn(s) 


( 36 ) 


The  required  factorisations  may  be  performed  approximately  if  we  assume  that 
AiCji  <<  1 .  By  using  the  sum  split*^ 

sgn(s)  «=  1  +  —  £n  s  -  —  Cn  s  ,  (37) 

where  the  dimensionless  variable  s  is  chosen  to  be  s  =  s/|si  and  £n+s  and  £n_s  have 
branch  cuts  in  the  low'er  and  upper  half-planes  respectively  and  are  real  when  s  is  real 
and  positive,  it  follow-s  that* 

1  +  i  Ac  .  sgn(  s )  -  [d.(l+iAc.(e.  +  —  £n  s )  )  j  •  d  .  1  +  i  Ac  .  ( 1  -e  .  -  —  Cn  s  ;  1  '  ■+ 

0  [  i  Ac .  1 2]  ;  i  -  1  ,  2  (38.) 

(dj,  d2.  ej  and  e2  are  arbitrary  constants.)  The  quantity  in  the  first  (second)  set  of 

squared  brackets  is  regular  as  s  ->  a  in  the  half-plane  Im  s  >  -ko  (Im  s  <  k^  cos''0). 
We  conclude  that 

f+(s)  and  g+(s)  -*  constant  as  s  -»  ^  in  Im  s  >  -k^  (39a) 


f_(s)  and  g_(s) 


constant  as  s 


or  in  I  m  s  <  k  cos  0 
o  o 


(39b) 


From  (29b)  and  (29d)  it  now  follow's  that 

S  (s)  and  T  (s)  0  as  s  ^  ^  in  Ims  >  -k1  .  (40) 

+  +  o 


In  accordance  with  our  earlier  remarks  we  see  from  (30a)  that  for  J(s)  to  vanish 
everywhere  on  the  complex  s-plane  it  is  sufficient  to  assume  that  either  p>- 2  or  i  >-l  . 
Similarly,  for  K(s)  to  vanish  identically  it  is  sufficient  that  either  p  >  0  or  <t>-1  These 
conditions  are  satisfied  if,  for  example,  we  require  that  Vy(xi>')  possesses  at  most  an 
integrable  singularity  at  the  point  x  =  y  =  0- 


For  the  moment  we  assume  a  unique  solution  to  the  problem  so  that  J(s)  and  K(s) 
are  zero  for  all  complex  s.  The  scattered  field  is  then  given  by 


oo+  i  c 


1 

— >  ( s ) | y | - i sx  r  T. 

(s) 

y  S  .  ( s  '  I 

4  T  \ 

1.  e  k(s) 

g-(s) 

i  y  i  f  _  ( s )  J 

-o>+  i  c 


where  we  have  assumed  that  f_(s)  and  g_(s)  are  non-zero  in  the  half-plane  Im  s  < 
k^cos^.  With  the  known  asymptotic  dependence  of  the  functions  S_(s),  T_(s).  f_(s)  and 
g_(s),  it  is  a  straightforward  matter^  to  deduce  that  p  ^  1  and  r  ^  0  for  this  solution. 
This  behaviour  confirms  that  the  set  of  sufficient  edge  conditions  derived  above  is  indeed 
satisfied.  Also*  it  ensures  that  the  component  of  the  total  field  parallel  to  the  wedge  (ie 
v?0)  is  finite  at  the  edge,  a  necessary  condition  on  the  solution  since  it  is  a  consequence 
of  the  requirement  that  the  electromagnetic  energy  in  any  finite  domain  be  finite^. 


*lf,  for  example,  we  choose  to  define  s  =  s/k£  then  the  terms  omitted  from  the  right 
side  of  (38)  are  small  only  for  tsi®k^ 
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From  (41)  we  see  that  factorisation*  of  f(s)  and  g(s)  are  required  for  a  determination 
of  v.  Since  J(s)  and  K(s)  are  now  known  to  vanish  identically,  (30)  may  be  used  to 
reuTite  (24a)  and  (24b)  as 


r+(s) r_(s) 


i  +  ia  (_L_  -  j,s) 

b  h(s)  ’  f_(s) 


2->  (s) 


(s-k  cos#  ) 
o  o 


f'(k  cosfl  ) 

i  ( s -k  cos 0  )  - — —2 - — 

o  o  f  (k  costf  ) 
_ _  +  o  o 

f '  (  k  c  o  s  0  ) 

(s-k  cos$  )  fX-  ° - 

o  o  i  v  k  c  o  s  H  ) 

+  o  o 


(42a) 


g+(s)g_(s) 


1  +  iA  - 
a 


[k?-K2] 

■) 3  ( s ) 


■(«) 


*  (a‘bh<!>  rui  * 


g'(s)  [2k+scos0  +  k  cos^fl  1 
,  l  o  o  oJ 


k+s  cos s-k  co s ^ 
oJ  l  o  o 


fk+k  cos^e  1 f s-k  cos«  1  g'(k  cos^  ) 
l  o  oil  o  oJ  °+  o  o 

[2k  +  SCOS0  4-  k  cos  0  |  2  (k  cosfl  ) 

l  o  o  0  J  °+  0  0 

(k+k  cos2^  ]  [s-k  co  s  $  ]  g’ (k  cos*  ) 
l  o  oJ  l  o  oJ  o  o^ 

[k+scosfl  ]  g  (k  cos*  ) 

l  oJ  0  +  o  o 


where 

bk2  -  ak2 

k  =  ~a-b)k  °  <4-4) 

o 

To  proceed  further,  it  appears  that  complicated  'transcendental'  factorisations  must  be 
found:  ie  the  functions  to  be  factorised  depend  upon  the  results  of  the  factorisations!  In 

fact  the  situation  is  not  so  serious.  If  we  restrict  the  analysis  to  those  thin-wedge 
(A  «  I)  problems  where  the  electromagnetic  parameters  are  such  that 


A  «  1 

b 


a|—  «  1 

a 


b  rr 

A  -  - T 

3  k2 
o 


«  1  (44) 
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then  approximate  factorisations  are  possible.  When  these  conditions  are  satisfied  it  is 
meaningful  to  look  for  solutions  f±  and  gt  of  the  form* 


ft(s) 


i  A 


(a-b) 

b 


(45a 


S,(s) 


4 


i  A 


(a-b) 

a 


0(A2) 


i  + 


iA^(s) 


g4>  +  0(A 


!>) 


(45b) 


where  f]  2  and  £1  2  3  4  are  arbitrary  finite  constants.  This  leads  to  two  simplifications. 
First,  the  functions  to  be  factorised  may  be  expanded  in  powers  of  the  small  parameters 
(44).  We  choose  to  neglect  those  terms  that  are  at  least  of  order  O(A^).  Second,  the 
factorisations  are  reduced  to  the  easier  problem  of  performing  sum  decompositions,  with 
the  errors  again  of  order  O(A^).  The  sum  splits  to  be  found  are 


F(s)  -  F+(s)  +  F_(s) 


S _  _  2«y  (s) 

( s )  ( s  -  k  cos^'  ) 

o  o 


( 46a ) 


C^s)  -  C+(s)  +  C  ( s ) 


(a-b)  .  2y (s) 


a  (s  -  k  cos^  ) 

o  o 


4 


hu2\\ 

s  2 

i,3(s)  ■>(s)(s  '  k0C0SVJ 

l  1  oJ 

cos^ 


j^s  cos^  +  kj 


(a-b) 


(s) 


7  (s) 


1 


(46b) 


Thus,  when  the  functions  F  +  and  G*  have  been  determined,  the  factorisations  (42)  will  be 
known,  correct  to  first  order  in  the  parameters  (44). 

The  decomposition  (46a)  is  carried  out  in  the  usual  manner.  We  find 


*These  are  the  most  general  forms  consistent  with  the  known  behaviour  (39)  of  the 
functions  as  isi  ®  in  the  appropriate  half-planes. 

Note  also  that  provided  we  choose  the  parameters  (44)  small  enough  we  can  validate  our 
earlier  assumptions  on  the  non-vanishing  of  f ♦  and  gt  in  the  appropriate  half-planes. 
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Ft  (s) 


*(s-k  cos'?  ) 
o  o 


2k  #  s  i  nf?  + 
o  o  o 


[2ko-s(s*kocos«o)] 


-1 

t  an 


s 


I 


(47) 


when  is i  <  | kQ i . 

A 

Although  it  appears  that  the  function  G(s)  is  singular  at  the  point 


k 

s  -  - 

cos  0 

o 

this  is  not  so  for  either  set  of  boundary  conditions  of  physical  interest.  When  a  =  1  and 
b  =  f  .  (  ,  (H-polarised  incident  plane  wave),  k  vanishes  and  the  singularity  is  removable. 
Also,  when  a  =  b  =  1  (E-polarised  incident  w'ave),  |ki  becomes  infinite  and  the  final 
term  on  the  right  side  of  (46b)  is  absent.  The  result  of  the  sum  split,  again  valid  for 
i  s  i  <  i  kQ  i .  is  as  follows: 


C  is  'i  - 


( a  -  b )  i 


cos'5 


s-k  cos'5  s  cos'5  +k 
o  o  o 


(k2-s2l 
t  o 


-1 

t  an 


5  I 


a 


2  f  k  0  sin'5 
o  o 


r ( s - k  cos'? 
o 


o 

~ 

o 


2 


r  s-k  cos'5 
o  o 


cos'? 

o 

scosfl  +k 


k2-S2 

o 


K-1] 


-1 

t  an 


s 


4 


T 

2 


b 

a 


20 

o 

k  sintf  (s-k  cos'?  ) 

r\  A  A  A 


1 


(48) 


(F  (s)  -f  F  (k  cos  0  )  +  (s  -  k  cosO  )  F '  ( F;  c  o  s  ^  )  )  1  )  (50 
+  o  o  o  o-fo  oJJ. 

Again,  expansions  in  the  small  parameters  (44)  have  been  used  to  arrive  at  this  form 
This  result  gives  exact  expressions  for  the  first  two  terms*  of  an  expansion  of  ■„  in 

powers  of  the  wedge  angle  A.  In  particular,  note  that  to  this  order  all  dependence  on 

the  arbitrary  constants  fj  i  and  £l  2  3  4  disappears. 

5  DERIVATION  OF  EDGE  DIFFRACTION  COEFFICIENTS 

Contour  integrals  of  the  type  (49)  arise  frequently  in  the  treatment  of  problems  to 

which  the  Wiener-Hopf  technique  may  be  applied.  Noble ^  has  given  a  useful  account  of 
the  analysis  of  this  integral. 

At  this  stage  it  is  convenient  to  let  tend  to  zero,  so  that  the  integration  contoui 
lies  along  the  real  axis  of  the  s-plane,  with  indentations  as  shown  in  Figure  2a  It  is 
standard  procedure  to  revert  to  polar  coordinates  and  to  make  the  substitution 

s  =  -kocos0.  Under  this  transformation  the  s-plane  is  mapped  into  an  infinite  number  of 
parallel  strips  of  the  0-plane  Choosing  the  segment  0  ^  Re£  ^  t  the  integration  contour 
becomes  the  contour  C  of  Figure  2b.  The  integral  becomes 

ik  rcos(0*'M 

v(r,P)  -  k  f  w(-k  cos0)  e  °  sin0  d£  ,  (5U 

o  J  o 

C 

w'here  the  upper  (lower)  sign  applies  for  0  ^  0  £  r  (-*  ^  0  4.  0). 


*It  is  easily  verified  that  the  inclusion  of  terms  of  order  O(A^)  in  the  abe’s  (13)  leads  to 
corrections  to  the  scattered  field  which  are  at  least  of  order  0(A*^) 
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Figure  2(a).  Integration  contour  in  the  s-plane. 


h 


Figure  2(b). 


Integration  contour  C  and  steepest 
descents  path  D( J 0 ' )  in  the  F-plane. 


An  approximate  evaluation  of  (51)  for  large  values  of  the  parameter  kQr  may  be 
obtained  by  the  method  of  steepest  descents.  The  saddle  point  of  the  integrand  is  at  Ci  - 
0$  =  \0\  and  the  path  of  steepest  descent  is  D(|fM).  (See  Figure  2b.)  In  the 

deformation  of  the  path  C  into  the  new  contour  D(  i  l ) ,  the  pole  at  0  -  r  -  is 

captured  when  1 0 1  >  *  -  For  0  ^  6  4.  the  threshold  occurs  at  0  =  t  -  oQ  and 

corresponds  to  the  reflection  boundary  of  the  half-plane  y  =  0,  x  <  0.  Thus  the 

contribution  of  the  pole  is  associated  with  the  reflected  wave.  Similarly,  when  -t  ^  0  £ 
0,  the  pole  term  is  included  only  in  the  shadow  sector,  and  gives  rise  to  the  field 
transmitted  through  the  dielectric.  The  integral  along  the  steepest  descent  path  D(  i  e  i )  is 
interpreted  as  the  diffracted  field  vd(r*^)-  Thus 


i  k  rcos (0- I  0 \ ) 


v(r ,0) 

=  V  ,  ( r , P )  +  mk  f  w(-k  cos0)  e  °  sin0  d0 

d  o  J  o 

CP 

(52 

with 

vd(r,P) 

ik  rcos (0- l 0 \ ) 

=  k  1  w ( “k  co s 0 )  e  s  i  n0  d /; 

o  J  0 

D(lPi) 

(5? 

and 

o  o 

1  1 

*  * 

V  A 

p 

11 

£ 

(54) 

The  contour  Cp  encloses  the  pole  at  0  =  r  -  (and  no  other  singularities)  and  is 

traversed  in  the  positive  (anti-clockwise)  sense. 

Since  we  are  primarily  concerned  with  a  study  of  diffraction  w’e  shall  consider  v'd 
only.  An  asymptotic  expansion  of  the  type  suggested  by  Jones^  may  be  used  to  evaluate 
v'd  for  kQr  >>  1.  However,  for  the  purposes  of  extracting  the  diffraction  coefficient,  it  is 
sufficient  to  determine  the  diffracted  far  field,  which  is  obtained  from  the  leading  term  of 
such  an  expansion.  Upon  inspection  of  (50)  it  is  apparent  that  this  procedure  breaks 
down  in  the  vicinity  of  the  reflection  and  shadow'  boundaries;  ie  when  the  pole  of  the 
integrand  at  0  =  r  -  is  close  to  the  saddle  point  at  0  =  1 0  i . 

Thus,  by  the  method  of  steepest  descents,  it  is  readily  shown  that 


k  r  -> 

l  ,(r,P)  — - >  /2x  k  isinfM  w(-k  cosf1, 

d  o  o 


i (k  r  -  t /4 ) 
o 
e 


7T. 


O 


r 


s  s 


provided  that  cosf1  +  c. 

As  expected,  in  the  far  field  v’d  appears  as  a  cylindrical  wave  originating  at  the 
diffracting  edge.  The  diffraction  coefficient  )  is  defined  by  the  relation 

k  r  ->  -  ikor 

v-d(r,fl)  _2 - >  D^,  0o) 

o 


From  (50),  (55)  and  (56),  then,  we  deduce  that 
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( a -b)  e 

D )  «  -A  - - -  — 

o  a 


c osO+cosfl 


[k-k^cosftcosfl  ^ 


[1  -  iA(C  (-k^cosfl)  +  C^Ck^cosfl^) ) ] 


+  iAfk+k  cos^  l(cos^cosfl  )  C’(k  cos 0  ) 
[  o  oJ  o+oo 


-  5  s  i  s  i  1 
b  o  l 


iA(a-b) 


A 

(F-<'ko 


cos/5)  + 


+  F  (k  cosP  )  -  k  (cosfl+cosfl  )  F'(k  cos^  )]])  (57) 

+  0  o  o  o  +  o  oJJJ 

for  0  £  f'Q  <  r-A  and  cosi*  +  cos90  0,  provided  that  the  strong  inequalities  (44)  are 

satisfied. 

It  is  clear  that  Cj  is  invariant  under  the  simultaneous  replacements  0  -*  -P  and 
"o  “*  -"o  This  *s  consistent  with  the  mirror  symmetry  of  the  diffracting  object  in  the 
plane  v  =  0  and  so  we  deduce  that  the  expression  (57)  obtained  for  the  diffraction 
coefficient  is  valid  also  for  -(r-A)  <  0o  ^  0. 

From  the  reciprocity  of  Maxwell’s  equations  it  follows  that  the  diffraction  coefficient 
should  be  unaltered  by  the  interchange  of  incident  angle  P0  and  diffraction  angle  0 . 
Indeed,  this  invariance  should  apply  to  each  term  of  an  expansion  of  the  diffraction 
coefficient  in  powers  of  wedge  angle  A.  It  remains  to  rewrite  (57)  in  a  form  which 
makes  this  symmetry  evident.  By  using  the  explicit  forms  (47)  and  (4S)  for  the  functions 
F  ♦  ( s )  and  G1**^(s),  it  is  found  that 


D;  'G'V 


cos/'+cos/1 


,  . .  (k  cosPcosO  -k) 

( a -b^  o  o 


i  .  -  (a~b)  frffl  o  1  +  rU  ! 

7T  a  *  o  to’ 


(a-b)  .  . 

+  A  — ; sin/1  sinfl 

b  o 


1  -  - 


A  (a-b) 


~  (t(P,0o)  +  t(PG  ,0)) 


\o 

,  0  1 

(k+k0C°S2^l ) 

20  si n0 

1  +  1  *  ( 

L  1 

2  J 

2  1 

k  sin 

o  1  1 

(cosfi  +cosfJ  ) 

k  0  sinO  cosf* 
o  1  \  2 

k-k^costf  ^ cosf^ 


*(*]  ,o2) 


icotpi(k+kJ 

k  sin 
o  I 


0 .  \  1  +  s i n  ^  .+cosff  1  cos^-1 

1  +  — _ _ ! _ 1  2J 

s\nO  (cos(?  +co sO  ) 


h  is  easily  checked  that,  once  the  differences  in  notation  have  been  reconciled  the 
SiT»?  Sri"  'h°“  l,n“r  'n  W'dE'  an*,e  A)  '«Prod“'  “*  '«“»•  oblained  t, 

6  DISCUSSION  OF  RESULTS 

lK  r2£*  *C/'n  b>  summarising  our  results.  When  the  incident  plane  wave  is  H-eolarised 
the  GTD  diffraction  coefficient  is  *■ 


D  0  ) 

o 


2 

cos^+cosf' 

0 


L 

-  ^  A  - - -  cos^cos^ 


A  (n  -11  r_H  h 

-  — r—  r  )  +  rH<*  p) 

n‘  1  °  ° 


A i  n  - 1 )  s  i  n^s  i  nf1 


[l  -  *  (n2-l  )(<t(P,eQ)  +  4(f>o,f>))]! 


where 


H  i  2f)  tan^  cosf< 

r  -  cot 2 ft  1  -  _I_ — 1 _ i  .  *  coxe 

1  1  j  (cosA j+cos^^)  rlCOlf] 


4(^1  ^2)  is  given  by  (60)  and  n  =  kj/k0  is  the  refractive  index  of  the  dielectric  relative 
to  the  surrounding  medium  (taken  to  be  air). 

For  an  E-poIarised  incident  wave  the  diffraction  coefficient  is 
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D  O.'’  )  - 

o 


A(n2-1) 


|cOS^+COSflo j ' 


1  _  A  (n2-l)(rE(^,(?  )  +  rE(^  (?))] 

It  O  O  J 


(63) 


with 


rE<VV 


cosec  0 


20  ^sinO  ^ 

1  +  - - - - — -  -  0.cot0t 

(cosO  +cos^)  1  1 


(64) 


In  both  cases  the  diffracted  far— field  is  given  by  (56).  These  expressions  are  valid  only 
when  the  order  O(A^)  corrections  are  small  in  comparison  with  the  leading  contributions. 
This  requirement  places  a  number  of  restrictions  on  the  permissible  values  of  the  various 
angles  and  parameters: 

(i)  -Assuming  that  n^  >  1,  (61)  and  (63)  are  useful  only  when  (n^-1)  A  <<  1 . 

(ii)  We  have  already  noted  the  breakdown  of  the  analysis  in  the  vicinity  of  the 
reflection  and  shadow  boundaries,  ie  when  Oqca  *  ♦  0. 

(iii)  We  observe  that  the  functions  and  $  are  all  singular  as  ♦  r, 

implying  that  the  formulae  (61)  and  (63)  are  invalid  when  either  the  incident  or 
diffracted  wave  approaches  one  of  the  faces  of  the  thin  wedge. 


It  is  worth  pointing  out  that  although  equations  (61  )  and  (63)  show  the  diffraction 
coefficients  as  expansions  in  powers  of 


n 


A 


2 

and  t  n* -  1  }  A 


there  is  no  guarantee  that  such  a  pattern  will  persist  to  higher  orders  The  missing  terms 
are  considered  negligible  only  because  of  the  assumption  that  A  <<  1 . 


The  results  (61)  and  (63)  remain  valid  for  a  lossy  dielectric  provided  that  the 
condition  i  n2-l  i  A  <<  1  is  not  violated  Therefore  the  results  of  this  analysis  cannot  be 
applied  to  a  studv  of  scanning  from  a  thin  wedge  with  high  conductivity  In  particular, 
a  reconciliation  with  the  known  dufraclion  coefficients  of  a  perfectly  conducting  wedge  is 
not  possible 


*Of  course,  a  similar  Wiener-Ffopf  type  approach  with  suitable  approximate  boundary 
conditions  may  be  used  in  a  study  of  scattering  by  a  thin  perfectly  conducting  wedge 
The  essential  difference  is  that  in  this  case  the  diffraction  coefficients  remain  non-zero  as 
A  becomes  vanishingly  small,  since  then  the  diffracting  object  is  a  perfectly  conducting 
ha  If -plane. 
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Kaminctzkv  and  Keller^  have  used  a  perturbation  approach  to  deduce  diffuct.  n 
coefficients  for  near!)  transparent  (diaphanous)  wedges  Specific  alls,  expressions  hue.o  m 
the  small  parameter  f,  which  is  a  measure  of  *he  difference  of  the  refractive  mdt  >  of  the 
dielectric  from  that  of  the  surrounding  medium,  have  been  obtained  F  urthermore  i:  was 
verified  in  Ref  2  that  both  this  type  of  analysis  and  a  perturbative  studs  of  thin  wedres 
give  identical  leading  terms,  le  terms  of  order  0(eA),  for  diffraction  coefficients  of  wedges 
which  are  both  thin  and  almost  transparent  It  is  a  simple  task  to  use  the  results  of  the 
present  work  to  extend  this  check  to  include  a  comparison  of  contributions  of  ordv: 
O(rA^).  In  fact,  it  is  found  that  in  both  types  of  expansion  there  are  no  such  terms 

It  has  already  been  remarked  that  the  leading  contributions  to  the  thin  wedge 

dielectric-edge  diffraction  coefficients  (61  )  and  (63)  are  in  accord  with  the  earlier  work  c  t 
Kaminctzkv  and  Keller  The  effect  of  including  the  second  order  terms  is  illustrated  in 
Figure  3,  where  the  magnitudes  of  the  diffraction  coefficients  for  a  normally  incident  plane 
wave  are  plotted  as  functions  of  wedge  angle.  The  corrections  are  more  significant  for 
the  case  of  E-polarisation,  where  they  are  of  opposite  sign  to  the  leading  contributions 

In  Figure  4  the  behaviour  of  the  diffraction  coefficients  is  shv*n  as  a  function  of 

observation  angle  for  various  angles  of  incidence 

It  appears  that  there  is  no  available  expennunt.il  data  on  the  scatter  ir.gr  of 
electromagnetic  radiation  by  penetrable  wedges  Yeo.  W.il!  and  Rates^^  have  embarked 
on  a  programme  of  work  aimed  at  developing  an  alg  :.:hm  foi  the  computation  of 

dielectric-edge  diffraction  coefficients  The  results  obtained  To  small  wedge  angles  uvrg 
this  combined  analytical-computational  approach  have  been  compared  with  those  of  Ref  2 
and  show  good  agreement 

7  CONCLUSIONS 

A  detailed  analysis  of  plane  wave  diffraction  bv  a  thin  dielectric  wedge  has  hr 
given.  Particular  attention  has  been  given  to  the  role  of  edge  conditions  and  then 
relevance  to  the  question  of  uniqueness  of  the  solution  The  solution  for  the  scatter  J 
field  has  been  obtained  in  the  form  of  a  contour  integral  The  fields  reflected  and 

transmitted  by  the  wedge  are  associated  with  the  pole  terms  In  the  fai  field  region  the 
remainder  of  the  solution  appears  to  emanate  from  the  edge  of  the  wedge  and  is 
interpreted  as  being  due  to  diffraction  by  the  edge  The  corresponding  GTI)  edge 
diffraction  coefficients  are  obtained  as  expansions  in  the  small  wedge  angle  A.  but  are 
valid  only  when  in^-1  i  A  <<  1  The  results  comply  with  the  requirements  of  reciprocr.y 
and  are  exact  to  order  O(A^) 
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